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Abstract. We consider 1+3 dimensional maximally symmetric Minkowski brane 
embedded in a 1 + 4 dimensional maximally symmetric Minkowski background. The 
resulting 1 + 3 dimensional effective field theory is of DBI (Dirac-Born-Infeld) Galileon 
type. We use this model to study the late time acceleration of the universe. We study 
the deviation of the model from the concordance ACDM behaviour. Finally we put 
constraints on the model parameters using various observational data. 
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1 Introduction 

Providing a satisfactory explaination of the late time acceleration of the universe is 
one of the most challenging tasks for cosmologists and particle physicists at present 
[1]. The most studied approach invloves adding an exotic form of energy with negative 
pressure (termed as dark energy) in the energy budget of the universe. The cosmo- 
logical constant ( with w = - = — 1) is the simplest candidate for such a component. 
Although the concordance ACDM model is allowed by all current cosmological ob- 
servations, it is also plagued by serious issues like fine-tuning and cosmic coincidence 
problems. Scalar field model [2] is another example of dark energy where the equation 
of state w evolves with time which in turn helps to solve the cosmic coincidence by 
tracker- type evolution. This also restricts the form of the potential for the scalar 
fields. 

Another interesting approach to explain the late time acceleration of the universe 
is to modify the gravity at large scale ( infra-red modification of gravity). DGP brane- 
world model [3] is one such example where the gravity is altered at large scales due 
to the slow leakage of gravitons from our observable universe ( modelled as an three 
brane) in to the higher dimensional bulk. The resulting Hubble equation can lead to 
late time acceleration in the observable universe. DGP model provides an interesting 
setup to modify gravity where one can construct a new four dimensional effective field 
theory which contains nontrivial symmetry properties. These symmetries are due to 
the combination of five dimensional Poincare invariance and the brane paramatrization 
invariance. It has been shown that in this effective theory, there is a single scalar field 
7r which represents the position of the three-brane ( our observable universe) in the 
higher dimensional bulk. The effective action contains a cubic self-interaction term of 
the kind (97r) 2 D7r together with normal canonical kinetic energy term. This cubic self 
interaction term has the property that it leads to second order equation of motion of 
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7r. This term is also invariant under the Galilean Transformation: 



tv — > tv + a 



df.Tv ->■ d^n + b 



(1.1) 



where a and 6 M are constants. Due to the form of equation (1.1) vr is often referred as 
Galileon field. Later, a four dimensional theory for the field tv with Galileon symmetry 
was proposed that contains five terms which also had the inersting property that depsite 
the presence of higher derivatives in the action, the equations of motions are second 
order [4, 5] . This causes the theory to be free of any ghosts. Cosmology with such an 
action has been widely studied in recent times [6]. 

On the other hand, Dirac-Born-Infeld (DBI) action contains the lowest order 
dynamics of a brane embedded in a higher dimensional spacetime. This gives an 
interesting setup to study inflation [7] and late time acceleration [8, 9] of the universe. 
If the universe is indeed described in a brane world scenario, then they should share 
the symmetries in the DBI action. Galileon terms can be thought of as a subset of 
all the higher derivative terms usually expected to be present in any effective field 
theory of the brane and they will be suppressed by some cut of scale (of different 
powers). In a recent work, Rham and Tolley [10] have constructed a general class 
of effective field theory which reproduces Galileon as well as the familiar DBI action 
under different limits. This has been extended by Goon et al. [11] where they have 
constructed a general class of effective field theory assuming that a 3-brane is moving 
in the higher dimensional bulk. Galileon theory arises as a special case of this setup. 
Similar generalization of Galileon field to maximally symmetric spacetime using a de- 
Sitter slicing has also been done in [12]. 

In this work, we study the late time acceleration of the universe in such a setup. 
We assume the simplest case where a maximally symmetric brane ( Minkowski Brane) 
is embedded in a maximally symmetric bulk (Minkowski bulk). We consider those three 
terms in the action which under small field limit generate the standard Galileon terms 
that arise in the DGP model under decoupling limit. We compare the cosmological 
behaviour in this setup to that of the concordance ACDM universe. We also constrain 
our model parameters using currently available observational data. 

We start by giving a general introduction on DBI-Galileon model in section 2. 
In section 3, we discussed the late time cosmology in this model. We explore the 
constraints on the model parameters using observational data in section 4. Section 5 
deals the conclusions. 



We begin with a setup where the foliation is Gaussian normal with respect to the bulk 
metric Gab and the extrinsic curvature on each foliation is proportional to the induced 
metric. This allows us to write the metric in the following form [11]: 



2 DBI-Galieleon 



G AB dX A dX B = f(Tv) 2 g^(x)dx^dx u + dn 2 



(2.1) 
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where ir represents the Gaussian normal transverse coordinate, and is the metric 
on the brane. We have choosen a gauge used in [11] 

X a (x)\ a=0 ... 3 = x» ] X 4 = tt(x). (2.2) 
Then the tangent vectors are given by 

where, 7r ;A1 denotes the covariant derivative of the field compatible to the metric g^ v . 
The normal vectors are given by 

^={7^=7 ^ 



; 1 1 (2-5) 

1 + 4(^7T ;/[1 7I>) 



where we have defined 7 as: 



The induced metric on the 1 + 3 dimensional hypersurface is defined as 

V = G AB e^. (2.6) 
In the gauge given in equation(2.2), the induced metric is 

hp, = fi^fg^u + n.^ii.v. (2.7) 

Let us now consider that a Minkowski brane (M 4 ) is embeddd in a Minkowski bulk 
(M5). We choose cartesian co-ordinates (x M , ir) on M5. Then, 7r = constant will give 
the foliation of M5 by M4 and the five dimensional line element can be expressed as 

ds 2 = n^dx^dx" + dir 2 (2.8) 

This is equivalent to putting 

/(7T) = 1,^=7^. (2.9) 

For this metric following Galileon terms are possible for which Galileon symmetry is 
preserved and equation of motion for the field ir is second order [10, 11]: 



Ci = 7T 

£2 = -y/l + Wf 

£ 3 = -[n] 2 + 7 V] 

U = -7([n] 2 - [n 2 ]) - 2 7 3 ([tt 4 ] - [nM) (2.10) 
£5 = -7 2 ([n] 3 + 2[n 3 ] - 3[n][n 2 ]) - 7 4 (6[n][vr 4 ] - e^ 5 ] - 3([n] 2 - [[n 2 ])^ 3 ]) 
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We have used the notation II for the matrix of second derivative of the field i.e. 
11^ = 7r ;w „. [IP] has been used to denotes the trace of powers of the matrix. For 
example [II] = 7r. ; ^. [II 2 ] = II^IP^. All the indices has been raised with the metric 
gV. We have used the notation [ir n ] to denote the contraction of powers of II with 7r ;ft . 
[vr n ] = 7r ;/i . U^ n ~ 2 > .7i>. For example, [vr 2 ] = vr.^, [vr 3 ] = tt. ii W / tt. u . Here also all the 
indices has been raised with the metric g^ v '. 

These terms in equation (2.10) are the so called DBI Galileon terms [10, 11]. 
Applying the small field limit to these DBI Galileon terms, one gets back the original 
Galileon terms [10, 11]: 





7T 




c 2 = 


-it*: 


) 2 


c 3 = 




) 2 [n] 


u = 




) ! (M 2 


£5 = 




i 2 ([n] 3 



i[n] [n 2 ] +2[n 3 ]) (2.11) 

The first three terms in the standard Galileon case, are the simplest ones rep- 
resenting a potential term, a canonical kinetic term, and a self interaction term that 
can also arise in the decoupling limit of DGP model. Cosmology in FRW background 
universe with these terms has been studied by various author [6]. 

3 Late time Cosmological Evolution with DBI-Galileon 



In our present study, we consider the corresponding three terms in the DBI Galileon 
case (£±,£2 and £ 3 in eqn (2.10)) and study the cosmology in a FRW background. 
We also replace the ir term in the £\ by a more general potential function V(ir). The 
corresponding Einstein's equations are given by 

3H 2 = p m + V(tt) + -jfi== (3.1) 



VT 



7T 



and, 

2H + 3H 2 = V(tt) + c 2 VT^ - (3.2) 



C37I" 7T 
1 -7T 2 

We have set 8ttG = 1. Equation of motion of the field tt is given by 



7T 3H7T 3H7TTC 

^ (1-^)3/3 + C2 7r^f + C3 (T^y + v{7[) - = (3 - 3) 



-A- 



where V(n)^ n denotes the derivative of the potential with respect to the field. C2 and 
C3 are the corresponding coupling constants for £2 and £3. 
We define the following dimensionless variables: 



x = 7r, y 



-"0 



V IV 

7!^' ~~JT V 
vv^ 



C 2 n 
a= Jp' P = — ■ 



(3.4) 

(3.5) 
(3.6) 



(Vn) 2 ' 

Here H represents the Hubble constant at present. In terms of these dimensionless 
variables, we form the following autonomous system of equation: 



x 



3Xy 2 h(l - x 



2\3/2 



3ax(l — x 2 



a 



+ ZphxVl - x 2 



(3.7) 



y =T 



Ax 
~ 3h 



a 



or 



3h 2 



(3x 2 (3\y 2 h(l- 
3h{l-x 2 ) \ 



2 Y' 2 - 3ax(l - x 2 ) 



a 



+ 30hxVl - x 2 



(3.8) 



1 + 



a 



vT 



or 



(3x 2 



3h 2 



A 



3/i(l - x 2 
\ 2 x 



3Xy 2 h(l - x 



2\3/2 



3ax(l — x 1 



a 



3f3hxVT 



x' 



h 



(r 



(3.9) 
(3.10) 



Here, prime denotes derivative w.r.t ln(a). The equation of state parameter io„ 
for the field 7r is given by 



3h 2 y 2 (l - x 2 ) + a{l - x 2 f' 2 - (5x 2 hx 
a + 3h 2 y 2 \/\ — x 2 



(3.1i; 



We evolve the system from the decoupling era (a ~ 10~ 3 ) to the present day 
(a = 1). In the beginning we assume that there was negligible contribution from the n 
field and the universe was dominated only by matter. We also assume that the field 7r 
was initially frozen due to large Hubble friction. This is similar to the thawing class of 
models previously studied for both scalar fields [13, 14] and galileon fields [15]. This 
sets the initial condtions for h and 



hi = lOVlOOno, Zi«0. (3.12) 
The initial condtion for y is set in the following way: using eqn (3.1), one can write 

Vl = - ... " (3-13) 



3h 2 ^l 



x 2 
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Figure 1. Plot of w n vs ln(a) (left) and Cl n (right). In each plot, the three sets are for 
Xi = 0.2,0.5,0.9 from bottom to top. In each set, V(tt) = ir,TT 2 ,e w and l/ir 2 from top to 
bottom. f2 mo = 0.3 and a = 0.3, /3 = 0.1. 



This allows us to replace yi in terms of other quantities. On the other hand, 
Q n i is related to Q^o which in turn is equal to 1 — Q m o due to the flatness condition. 
Hence yi is not an independent quantity but is related to other model parameters like 
a, and £l m o- The initial value for A, Xi, is an interesting parameter. It controls the 
deviation from the ACDM behaviour. This is shown in Figure 1. It is evident from 
this figure that for smaller values of A, models with different potentials can hardly be 
distinguished from each other as well as from w = — 1. As one increases Xi, they start 
deviating from w — — 1 and deviate from each other as well. Moreover the equation of 
state Wyr, for linear potential has the largest deviation from w = —1. Note that this 
potential arises naturally in the DBI Galileon model. This result is also similar to the 
standard scalar field [14] as well standard Galileon field models [15]. 

One can also Taylor expand the scale factor of the universe around the present 
era (t = to) as follows: 

a(t) = a(t ) + a(t ) [H (t - t )] n , (3.14) 

n. 

n=l 

where, 

In = ^/^H n ). (3.15) 

It is now straightforward to show that —72 = q is the deceleration parameter. Similarly 
73 is related to the Statefinder r or jerk j and 74 is related to the snap s and so on. 
Using the equations (3.14) and (3.15), one can now construct the Statefinder Hierarchy 



- 6 - 



n m0 = 0.25, A, = 0.9 




1.0 



1.0 



n m0 = 0.3, Xi = 0.9 




0.2 



0.4 



0.6 



0.8 



1.0 



n m0 = 0.25, A, = 0.6 



n m0 = 0.3, Xi = 0.6 





1.0 



0.2 



0.4 



0.6 



0.8 



1.0 



Figure 2. Phase space diagram in S2 — plane for different potentials. V(ir) = r K,-K 1 ,e K 
and from bottom to top. a = 0.3 and /3 = 0.3 for all the plots. The vertical line represents 
the present day (z = 0). 



S n [16] as: 

5 2 = 72 + \tt m (3.16) 

53 = 7s (3.17) 

5 4 = 74 + yfi m (3.18) 

5 5 = 7B-3fi ro -yJ4 (3.19) 

33 34 
^6 = 76 + y^m + 3 4 fi^ + — fi^ and so on. (3.20) 

It is interesting to note that for ACDM, S n = 1 throughout the history of the Universe. 
Hence any deviation from S n — 1 represents models with non-ACDM behaviour. S n 's 
can also be used to study degeneracies between different dark energy models. We 
now use these Statefinder Hierarchies to study DBI-Galileon models with differenet 
potentials. 
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Figure 3. Phase space diagram in S3 — f2 m plane for different potentials. V(ir) = r K,-K 1 ,e K 
and 4y from bottom to top. a = 0.3 and /3 = 0.3 for all the plots. The vertical line represents 
the present day (z = 0). 



In Figure (2), we show the behaviour in the 62 — fl m plane for different potentials 
and for different values of Aj. As in the case of the equation of state, here too the 
deviation from the ACDM is small for smaller Aj, and large for bigger Aj. Also the 
deviation from the ACDM behaviour is highest for the linear potential. The degenera- 
cies between models with different potentials is higher for smaller Aj. Also for smaller 
values of f2 m n, deviations from the ACDM behaviour as well as deviations between 
different potentials are higher. In this figure, we keep both a and (3 fixed at 0.3. For 
other values, one can get similar behaviour. The dependence on a and is discussed 
below. 

Next we study the evolution for different potentials in the S3 — Q m phase-space 
in Figure (3). For this, we also fix a = 0.3 and /3 = 0.3. This figure is interesting. 
Here for Aj being 0.9 and 0.6, different potentials deviates more from each other at 
present(z = 0) for A; = 0.6. This is opposite to what happens in the case of S2 — Vt m 
phase space where the deviation at present is higher for Aj = 0.9. We have also 
demonstrated that all the potentials converge to ACDM for very small values of Aj 
(Aj ~ 0.1) for which all S n 's are exactly 1. This shows that there exists some Aj which 
optimizes the deviation in S3 for different potentials. 
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Figure 4. Phase space diagram in S3 — S2 plane for different potentials. V(tt) = 7^,ir 2 ,e n 
and \ from top to bottom at the left side of each figure, a = 0.3 and j3 = 0.3 for all the 
plots., The dots represents the present day [z = 0). 



In Figure (4), we show that phase space is S3 — S2 plane. Here also the models 
with different potentials deviate more for smaller Q m o and higher Aj. the deviation 
between models in this phase space is also higher than the two phase spaces. 

In Figure (5), we show the dependence on the two parameters a and /3. In this case 
we fix A, = 0.9 and f2 m o = 0.3. This is because, this choice maximised the deviations 
between different potentials as shown in Figure (4). Figure (5) shows that deviations 
between different potentials are larger for smaller values of both a and (3. For some 
choices of a and /3, there is complete degeneracy between some of the potentials. This is 
evident for the case a = 0.3 and (3 = 0.8 where the linear and the quadratic potentials 
are fully degeberate. 

From all these phase space diagrams, one can conclude that different potentials 
deviate maximum from each other for smaller Q m Q and larger Aj. But this can be 
adjusted by suitably choosing a and (3 to make them completely degenerate. 
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Figure 5. Phase space diagram in S3 — S2 plane for different potentials. V(tt) = n,ir' 2 ,e n 
and \ from top to bottom at the left side of each figure. \ = 0.9 and f2 m o = 0.25 for all 
the plots. The dots represents the present day (z = 0). 



4 Observational Constraints 

With this in mind, we now put constraints on these four parameters (Aj, Q m o, a and 
/3) by using recent observational data. 

To start with, we consider the Type la supernova observation which is one of the 
direct probes for the cosmological expansion. In this case, one measures the apparent 
luminosity of the supernova explosion from the photon flux received. Cosmologically 
the relevant quantity is the luminosity distance cLl(z) denned as: 

The distance modulus \x (which is an observable quantity) is related to the lumi- 
nosity distance as 

u = m-M = 51og-^ + 25, (4.2) 

Mpc 

where m and M are the apparent and absolute magnitudes of the Supernovae respec- 
tively. We consider the latest Union2.1 data compilation [17] consisting of 580 data 
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points for the observable /x. 

Next, we use the observational data on Hubble parameter as recently compiled 
by Moresco et al. [18] in the redshift range < z < 1.75. The sample contains 19 
observational data point for H(z) spanning almost 10 Gyr of cosmic evolution. These 
values are given in Table 1. It also contains the latest measurement of the Hubble 
constant H [19]. 
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Table 1. H(z) measurements (in units [kms 1 Mpc x ]) and their errors. 



Lastly, we consider the combined BAO/CMB constraints as recently derived by 
Giostri et al. [22]. 

We start defining the comoving sound horizon at the decoupling as: 

' (4.3) 



VZJo a 2 h{a)^l + {3Q b0 /4Q l0 )a 

where fi 7 o and Qbo are the photon and baryon density parameter at preset respectively. 
We have assumed the speed of light c = 1 in our calculations, z* is the redshift at 
decoupling and is given by the formula derived by Hu and Sugiyama[23]. Taking the 
data from WMAP7 [24], we put z* = 1091 exactly. We also fix the redshift of the drag 
epoch at z d m 1020. 

We next define the acoustic scale: 

U = vr — , (4.4) 
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where g?a(^*) = H J** dz'/h(z') is the comoving angular-diameter distance. The dila- 
tion scale is also defind as [25]: 



D v {z) := [d 2 A (z)H Q z/h(z)} 



1/3 



(4.5) 



The 6dF Galaxy Survey [26] and the WiggleZ team [27] have measured this quan- 
tity at 0.106, z = 0.44, z = 0.60 and z = 0.73. Percival et al. [28] also measured 
at z = 0.2 and z = 0.35. One can combine these masurements with the WMAP-7 
measurment of l a [24] and can obtain the combined meaurement of BAO / CMB for the 

dA( - z * S) — i This has been given in [22]. From this one can obtain [22]. 



quantity 



Dv(zBAo) 



where 



Xbao/cmb z 
( d A {z 



(4.6) 



X 



£V(0.106) 

<Ia{z*) 

ZV(0.2) " 
dA\z*) 
ZV(0-35) 

D v (0M) 

D v (0.Q) 
d A {z*) 
\ D v (0.73) 



and 



/ 0.48435 
-0.101383 
-0.164945 
-0.0305703 
-0.097874 
\ -0.106738 



-0.101383 

3.2882 
-2.45497 
-0.0787898 
-0.252254 
-0.2751 



-0.164945 
-2.45497 
9.55916 
-0.128187 
-0.410404 
-0.447574 



- 30.95 
17.55 

- 10.11 

- 8.44 
-6.69 
-5.45 



-0.0305703 -0.097874 -0.106738 \ 
-0.0787898 -0.252254 -0.2751 



(4.7) 



-0.128187 
2.78728 
-2.75632 
1.16437 



-0.410404 
-2.75632 
14.9245 
-7.32441 



-0.447574 
1.16437 
-7.32441 
14.5022 



/ 

(4.8) 

is the inverse covariance matrix. The correlation coefficients for the r s /Dy pair of 
measurements at z = (0.2,0.35), z = (0.44,0.6) and z = (0.6,0.73), respectively have 
been given by [27, 28]. 

The results are shown in Figure (6) and Figure (7) where we show the 2a confi- 
dence contours in the f2 m o ~ \ plane for different values of a and 0. Dependence on 
a is such that for smaller values of a, there is an upper bound on Aj which makes the 
cosmological evolution very close to the ACDM. But as one increases a, this upper 
bound shifts towards the higher values of Aj, thereby allowing large deviation from 
ACDM behaviour. This is true for all the potentials. The allowed deviation from the 
ACDM behaviour is highest for inverse- squared potential. 
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Figure 6. Allowed region in 2a cnfidence level in the m o — Aj plane. Top ones are for Linear 
potential while the bottom ones for quadratic potential. The left ones are for a = 0.1 (solid), 
a = 0.5(dashed) and a = l(dotted) keeping = 0.1. The right ones are for = 0.1 (solid), 
(3 = 0.5(dashed) and /3 = l(dotted) keeping a = 0.3. 



On the other hand, if one varies (3 keeping a fixed, there is also an upper bound 
on Aj. But as one increases the value of j3, this upper bound shifts towards the smaller 
value of Aj forcing the models to behave very similar to ACDM. This is opposite to 
the a dependence. One should note that a and /3 control the contribution from the 
£2 and £3 part of the total action. Hence bigger contribution from £2 term allows 
larger deviation from ACDM whereas bigger contribution from £ 3 restricts the model 
to behave more close to ACDM. 

5 Conclusion 

In this work, we investigate the late time evolution of the universe in a DBI-Galileon 
Model where a Minkowski Brane is embedded in a Minkowski bulk. To simplify our 
analysis, we only keep the first three terms in the total action which under weak-field 
limit reproduce the standard Galileon terms present in the decoupling limit of DGP 
model. We also keep a general potential term V(n) in the action instead of standard 
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Figure 7. Same as Figure (6) but the top ones are for exponential potential while the 
bottom ones are inverse squared potentials. 



linear term and study different choices for V(%). We assume the field n to be initially 
frozen due to large Hubble damping and that it behaves like a cosmological constant. 
However, with time the field slowly thaws out from the frozen state and starts deviating 
from w = —1. The deviation depends on the initial value of te slope of the potential, 
Aj. For smaller values of A,, the evolution for all the potentials remain close to ACDM 
throughout. As one increases the value of A;, the evolution start deviating from ACDM. 
We study the degeneracies for the different potentials using the statefinder hierarchies. 
Finally we constrain our model parameters using the recent observational data. We 
show that larger contribution from £ 2 P ar t results in larger deviation from ACDM 
behaviour whereas larger contribution contribution from £3 part restricts the models 
to behave more close to ACDM. This is true for all the potentials. 

Although we have not studied the complete action for the DBI-Galileon model, 
still this study gives some interesting observational consequences for first three terms 
of the full action. It will be worthwhile to study the observational consequences with 
full DBI-Galileon action and this will be our future aim. 
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